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"^^ ' Abstract 

,J^ ' We present an algorithm to compute the Weierstrass semigroup at a 

jrt , point P together with functions for eacli value in the semigroup, pro- 

ai vided P is the only branch at infinity of a singular plane model for the 

, , ' curve. As a byproduct, the method also provides us with a basis for the 

, spaces C{mP) and the computation of the Feng-Rao distance for the cor- 

^— ( ' responding array of geometric Goppa codes. A general computation of 

^ ' the Feng-Rao distance is also obtained. Everything can be applied to the 

lO . decoding problem by using the majority scheme of Feng and Rao. 
•/i ' Key words - algebraic curves, singular plane models, approximate 

^~~] [ roots, Abhyankar-Moh theorem, Weierstrass semigroups, one-point Alge- 

^—^ . braic Geometry codes. Apery systems and Feng-Rao distance. 

^ . 1 Introduction 

j^ , One of the main problems in the theory of algebraic Geometry codes (AG codes 

^ ' in brief) is the expHcit computation of bases for the spaces C{G), G being a 

• • , divisor over an algebraic curve. There are some general methods to do this 

^<\ ' computation, like Coates or Brill-Noether algorithms. In this paper we study 

K> , an alternative procedure for the particular case G — niP, where P is a rational 

^ ' point of the curve, using the structure of the Weierstrass semigroup, the well- 

known theory of Abhyankar-Moh and the normalization of the curve. This will 
have also the advantage that it can be used for effective decoding of those codes. 
Thus, our aim is to review and develop links between coding and singularity 
theories. 

Our method is algorithmic and it assumes to have a singular plane model for 
the curve with only one branch P at infinity which is defined over the base field. 
This branch is nothing but the point P in the support of G. A first step uses the 
Abhyankar-Moh theorem to give a subsemigroup of the Weierstrass semigroup 
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r and functions achieving its elements, with the aid of the so-called algorithm 
of approximate roots. A second step computes the rest of the semigroup and 
corresponding functions by a triangulation algorithm starting from an integral 
basis for the integral closure of the algebra of the afSne part of the curve. 

In the case of curves over finite fields, a byproduct of our algorithm can 
be used for decoding the Algebraic Geometry codes C{m) with divisor of type 
G = mP by the method of Feng and Rao iQ. This method, based on a 
majority voting scheme, needs the knowledge of the Weierstrass semigroup at 
P and functions achieving its values, and it decodes errors up to half the so- 
called Feng- Rao distance 5FR.{rn) (an estimate for the minimum distance which 
is better than the Goppa distance). 

The integer SpRim) can be defined in terms of arithmetical relations among 
elements in the semigroup F. The precise value of dpRim) is known for few 
classes of semigroups as, for instance, for (many elements of) telescopic semi- 
groups (see ||l5[). Telescopic semigroups are complete intersection, and so a 
very special kind of symmetric semigroups (i.e. those such that c = 2g, c being 
the conductor and g the number of gaps). The Abhyankar-Moh semigroup is 
telescopic, but the Weierstrass semigroup can be quite general. 

The last part of the paper is devoted to the computation of the Feng-Rao 
distance for numerical semigroups S. Our computational method assumes that 
S is presented by means of its so-called Apery systems. Such systems are the 
natural way to describe semigroups when one is dealing with problems involving 



relations (see B). We show a formula (theorem 4.5) to compute SpRirn) for a 
general S. If S is symmetric, this formula is improved (theorem 4.6) for some 
elements with m > c. Also for symmetric semigroups, we show how the formula 

5FR{m) = min{r e5|r>TO + l- 2g] [@] 

holds for most values m E S with c > rn > 2c — 1 and, moreover, a formula for 
the minimum element rriQ e S such that the formula [@] holds for m > toq . An 
estimate for mo was given for telescopic semigroups in [|l5[ . 

Finally, we show how our algorithmic method to compute Weierstrass semi- 
groups r gives also, as a byproduct, an easy way to compute an Apery system 
for such semigroups and, hence, the Feng-Rao distance for their elements. 

2 Weierstrass semigroups and AG codes 

The following section outlines the connection between coding theory and Weier- 
strass semigroups, and it is abstracted from |1J] and |16). Consider a non- 
singular projective algebraic curve x defined over a finite field W such that x 
is irreducible over IF. One-point geometric Goppa codes are constructed as fol- 
lows. Take n + 1 different IF-rational points Pi, . . . ,Pn and P of the curve X: 



and take a positive integer m. Thus, one can consider the linear map 

evD : C{raP) — > F" 
/k^(/(Fi),...,/(P„)) 

and define the hnear code C{m) as the dual space of Im {evjj), whose length 
is obviously n. Denote by k{m) and d(m) the dimension over IF and the min- 
imum distance of the linear code C{m) respectively, where the integer d{m) is 
the minimum value of non-zero entries of a non-zero vector of C(rn). Goppa 
estimates for k(m) and d{m) are derived from the Riemann-Roch formula as 
follows (see ^^ for more details). If 2g — 2 < m < n, then 

k{m) = n — m + g — 1 
dim) > m + 2 - 2g ^ d* (m) 

where g is the genus of the curve x and d* (m) is the so-called Goppa designed 
minimum distance of C{m). Apart from the excellent asymptotical behaviour 
of their parameters (see pi| , for instance) , the main interest of these codes is 
that they can be decoded efficiently by the majority scheme of Feng and Rao 
algorithm, which will be briefly described below. 

Fix a rational function fi £ ]F(x) with only one pole at P of order i for 
those values of i for which it is possible, i.e. for the non-negative integers in the 
Weierstrass semigroup F = Tp of x at P. For a received word y = c + e, where 
c G C{m), one can consider the unidimensional and bidimensional syndromes 
given respectively by 

n n 

k=l fc=l 

Notice that the set {fi \ i < m, i e F} is actually a basis for C{mP) and hence 
one has 

C{m) = {y e F" I s,(y) = for i < m} 

Therefore one can calculate Si(y) from the received word y as 



St{y) = ^VkftiPk 



fc=i 

and thus Si{y) is called known for i < m. 

In fact, it is a known fact that if one had enough syndromes Sij{y) for 
i + j > m one could know the emitted word c, and all the involved syndromes 
can be computed by a majority voting (see [|lO|). 

This decoding algorithm corrects up to half the so-called Feng-Rao distance 



SpRifn) of C{m) (see definition 4.1 below). More precisely, the procedure cor- 

SpRi'lTT'') — 1 / • r n I IT 

rects up to errors, where m = nunjr G I | r > m\. In particu- 
lar, one has d{7n) > Spp{m') (see M] for a direct proof of this inequality). 



Notice that SpRim) only depends on the semigroup T. On the other hand, 
it is known that SFRim) > m + 1 — 2g = d* {m — 1) for aU m E T, the right 
hand side term being the Goppa estimate for the minimum distance. Thus, the 
method gives an improvement of the number of errors that one can correct for 
one-point AG codes. 

In practice, the main problem is computing F and the functions fi achieving 
the values of the semigroup F in order to carry out this decoding algorithm. On 
the other hand, a second basic problem is calculating the value of Spii{m) if one 
wants to know the number of errors that one can correct. The aim of this paper 
is to give an approach to the above two problems with certain conditions. 

3 Computing Weierstrass semigroups 

In this section we show how the Weierstrass semigroups of curves having a 
singular plane model with only one branch at infinity can be computed at the 
same time as functions achieving their values. Everything is based on the theory 
of Abhyankar and Moh about such plane models, and so this section can be 
regarded as an application of singularity theory to coding theory. 

3.1 Semigroups of plane curves with only one branch at 
infinity 

In the sequel, we will assume that IF is any perfect field, since the hypothesis 
on the finiteness of W will not be used. Nevertheless, we are always thinking 
of finite fields, because of the applications in coding theory. Let x be a non- 
singular projective algebraic curve defined over W such that x is irreducible over 
IF. Let T be a plane model for x with only one branch at infinity, i.e. such that 
there exist a birational morphism 



X 



T C p2 



and a line L C P^ defined over W such that L n T consists of only one geometric 
point P and x ha-s only one branch at P. Notice that, a fortiori, both P and 
the branch at P arc defined over the underlying field W, since T does. Thus 
there is only one point of x lying over P, which will be denoted by P. 

Set C ~ x\ {P} and Z = T \ {P}. One has the two following additive 
subsemigroups of IN: 

Tp = {-vpif) I / G 0^{C)} 
Sp - {-vpif) I / e OriZ)} 

The first one is nothing but the Weierstrass semigroup of x at P; this semi- 
group contains the second one but they are different unless the curve T is 
non-singular in the affine part. Notice that both IN \ Tp and IN \ 5'p are finite. 



In fact IM\rp has g elements, g being the genus of x, which are the Weierstrass 
gaps. In order to compute the cardinahty of IN \ Sp, we need the foUowing 
effective version of a standard fact in singularity theory. 

Proposition 3.1 Let A — 0^{C) and B = Ot{Z) be the respective affine 
coordinate W-algebras of the curves C and Z; then one has 

tt(rp \ Sp) = dimw{A/B) = Y, ^Qi^) 

Qez 

where Q ranges over all the closed points of Z and Sq{Z) — dimTp{Or ,q / Ot ,q) , 
Or,Q being the normalization of the ring Or,Q ■ 

Proof : 

The second equality follows from the fact that 

A/B^ Ot,q/Ot.q 

where Sing Z is the set of (closed) singular points of -E. In order to prove the 
first equality, we take an F-basis {hi, . . . , kg} of A/B, which can be computed 
with the aid of the integral basis algorithm [^ (see ||l^, |20) and also Q). 

Set B^ = B + Whi + . . . + Whi, for < i < s. We proceed by induction, 
so let < i < s and assume that we have found functions gi, . . . , gi which are 
linearly independent over IF such that 

B + Wgi + ... + Wg, = B' 

and 

-'^{9]) i Tp"^ Vj < i 

where V'p = SpV^ {-z>p(.9i), . . . , -v^{gj)\ Q Tp . 

Now look at /li+i . If — i>p(/ii+i) ^ Fp, then set g^+i — hi^i and go on. 
Otherwise, there exists f ^ B^ with w-p(/ii+i) — ifp{f) and such that 

-Vp{hi+i - /) < -Vj;{hi+i) 

Then we repeat the process with /i^+i — / replacing to /li+i . Since one has 
/ij+i ^ B*, it follows that in a finite number of steps we will be able to replace 
hi+i by 

gi+i = hi+i {mod B^) 

such that 

-Vpi9^+l) i r^ 

^ After |13|, the integral basis algorithm is efficient whenever the ramification of B\W is noj; 



wild. For our purpose, this procedure will be applied in the hypothesis of proposition 3.6 
what actually implies that the ramification is tame. 



At the end of the procedure s different elements in Tp \ Sp will be added, 
and hence {[(Fp \ Sp) > dim^{A/B). On the other hand, since A = i?' = 
B + IFgi + . . . + Wgs one has that any h £ A can be written in a (unique) 
way as h — g + Xigi + . . . + XsQs with g ^ B and A^ G F, the values v-p{gi) 
and v^{g) being pairwise different. Thus, one has either —v-p{h) e Sp or 
—ihp^h) = —Vp{gi) for a unique i. This proves the equality. 

D 

Now our aim is as follows: we intend to describe Tp by first computing Sp 
and the associated functions, and secondly by completing the semigroup up to 
Tp with the corresponding functions. According to the above proposition, this 
last thing can be done by means of the following procedure. 

Algorithm 3.2 (Triangulation algorithm) 

Input: Sp , Or{Z) and {hi, . . . , hg}. 

• Initialize S — Sp and B = Ox {Z) 

• For i = 1 , . . . , s do 

— Set gi =^ hi 

— While ^v-p{gi) G S do 

* Find f (z B such that ~Vp{gi — /) < ^^{gi) 

* Set gi ^ g.i ~ f 

— Set B = B + Wg, and S ^ S U {~Vp{g,)} 

• Next i 
Output: S and B. 

Remark 3.3 One can substitute the set S — S U {—iJ-pigi)} by the semigroup 
S — S — iJ-p{gi) ■ IN in the above algorithm, i.e. in each step one can cover more 
than one gap and also add the corresponding functions (products of existing data 
functions) to B. It may yield in general a faster algorithm, since one could stop 
just when I new values are added. This idea can be done effective after the 



results of the next section by using Apery systems (see remark 4-^^) 



3.2 Approximate roots 

The remaining part of the algorithm, that is, the description of the semi- 
group Sp and the construction of the corresponding functions, is known from 
the Abhyankar-Moh theorem and the so-called algorithm of approximate roots 



which, for the sake of completeness, will be restated below. Both were intro- 
duced and developed by Abhyankar and Moh in M. In particular, it will be 
shown that this way to compute the Weierstrass semigroup and the functions 
is effective and, moreover, the arithmetic properties of the involved semigroups 
will allow us to compute the Feng- Rao distance with the aid of some extra tech- 
niques, what will be done in section 4.3. First we need the notion of approximate 
root. 

Thus, let 5 be a ring, G € ^[Y] a monic polynomial of degree e and F G S'[Y] 
a monic polynomial of degree m with e\m. If we write m = ed, then G is called 
an approximate d-th. root of F if deg {F — G'^) < m — e = e{d — 1). The key 
result is that, provided d is a unit in the ring S, then there exists a unique 
approximate d-th root of F, which will be denoted app{d,F). You can see a 
constructive proof in ]^ , where it is shown in particular that the computation 
of approximate roots is very efficient. From now on, we will work with the 
coefficient ring S* = F[X]. 

Let the affine plane model T having only one point at infinity be given by 
the equation 

F = i^(X, Y) = Y™ -H ai(X) Y™-i + . . . + a„,(X) 

where m is actually the total degree of the polynomial F, and set n = degxF. 
In the sequel, we will assume moreover that 

(H) charW does not divide either degT or ep(T) 

With the above notations, one has m = deg T and n — deg T — ep(T), and thus 
(H) is equivalent to say that p = charTF does not divide either m or n. 

In this case, we can actually assume that p does not divide m = degT. In 
fact, if TO is a multiple oip but n is not, we choose k not divisible by p such that 
nk > TO,, and by doing a change of variables in the form X' = X -I- Y'^, Y' = Y 
we get a new affine curve which is isomorphic to the original oneg but whose 
degree is not divisible by p. This will be assumed in the sequel for simplicity, 
and thus we will be able to compute any approximate root of F. 

By the sake of economy one wants to avoid the computation of parametric 
equations for the singularity at infinity, and thus computing the semigroup 
Sp and corresponding functions directly from the equation of the curve. The 
parametrization would allow us to compute easily i>p-values, but this can be 
avoided by using resultants and approximate roots, as it is shown below. 

In order to state the algorithm of approximate roots from H , we first agree 
to set 

degxResY{G,H) = -cx) if ResyiG.H) = 

for any couple of polynomials G,H E IF[X, Y], and 

gcd {6o,Si,. . . ,Si) ^ gcd {6o,Si,. . . ,Sj) 



^In particular, Sp does not change. 



a 5o,Si, . . . , Sj are integers, j < i and Sj+i — Sj+2 — ■ . ■ — Si — — cxd. Thus, the 
algorithm works as follows with an input F as above (the case when Y divides 
F being trivial, we assume the contrary). 

Algorithm 3.4 (Approximate roots) 

• Set do = 0, Fq — X, So = di — m, Fi = Y and Si = degxResyiF, Fi) 

• For i from 2 do 

- di = gcd {di^i,Si-i) 

— If di — di-i then h — i — 2 and STOP else 

=1= F, ^app{d,,F) 

* St ^degxResY{F,Fi) 

• Next i 

Output: h, ((5o, ■ • ■ , 5h) and (Fq, . . . , Fh) 

Since the sequence {(ii}i>i obtained in the above algorithm is a decreas- 
ing one of positive integers, there exists a unique positive integer h such that 
di > . . . > d fi^i = dh+2 : and hence the algorithm stops. The first application 



of algorithm 3.4 is the following result, proved by Abhyankar in g], which pro- 
vides a criterion for a curve with an only (rational) point at infinity to have 
only one (rational) branch at this point (and to be absolutely irreducible, as a 
consequence) . 

Proposition 3.5 (Criterion for having only one branch at infinity) 

Let F he the equation of a plane model with an only point at infinity as above, 
and assume that char IF does not divide ni — deg F. Let h, di and Si the integers 
which are computed by the algorithm of approximate roots. Then the curve has 
an only (rational) branch at infinity if and only if dh+i ~ 1, Sidi > (^2^2 > . . . > 
Shdh and UiSi is in the semigroup generated by Sq,Si, . . . ,(5i_i for I < i < h, 
where Hi = di/di+i also for 1 < i < h. 

The second application of the algorithm of approximate roots is just the 
computation of Sp and the associated functions by means of the Abhyankar- 
Moh theorem, which provides us with a set of generators for Sp with nice 
arithmetic properties. The proof is refered to 0] or p^ . 

Proposition 3.6 (Abhyankar-Moh theorem) 

Let T be a plane model with an only branch at infinity, and assume that charW 
does not divide degT. Then there exist an integer h and a sequence of integers 
Sq, . . . ,Sfi G Sp generating Sp such that Sq — deg T and 



(I) dh+1 = 1 and Ui > 1 for 2 < i < h, where di = gcd (Sq, . . . ,Si-i) for 

1 < i < h + 1 and Ui = di/di-^-i for I < i < h. 

(II) UiSi is in the semigroup generated by 6o, . . . , 6i-i for 1 < i < ^ n. 

(III) UiSi > (5j+i forl<i<h-l. 



Remark 3.7 Notice that the restriction on the characteristic ofW is necessary, 
since the plane curve over IF2 given by the equation 

Y^+Y = X2(X*+X) 

has no affine singularity and an only point at infinity whose Weierstrass semi- 
group is generated by the elements {8,10,12,13}. This example gives a neg- 
ative answer to a question proposed by Pinkham in f^J^, since the sequences 
{8,12,10,13} or {12,8,10,13} actually satisfy (I), (II) and (III), but Jq = 8 
or (5o = 12, respectively, is not the degree of the curve. Thus, a general descrip- 
tion for T p = Sp in this way is still an open problem. 

Curiously, this example comes from coding theory (see uJ]), and thus the 
reason why the Abhyankar-Moh theorem fails is that both m and ep(T) are 
multiple of 2. 



Remark 3.8 Semigroups as in proposition p'.6| , as well as semigroups of values 
studied in positive characteristic by Angermilller ^in ^ or Campillo in /j^, are 
a particular case of telescopic semigroups, where only the properties (II) and 
dfi+i = 1 are required (see UW)- Being telescopic is equivalent to be free in the 
sense of M], what means that every m G Sp can be written in a unique (and 
effective) way in the form 

h 

m — 2^ ^i^i [*] 

i=0 

with Aq > and < A^ < n^ for 1 < i < h (see j7q// From property (II) one 
has that these semigroups are complete intersection (i.e. those such that the 
affine curve defined by them is a complete intersection one) and, in particular, 
they are symmetric (see /Ua/, /Uq/ and section 4.1 below for further details). 



The functions which are needed to achieve the pole orders given by the 
Abhyankar-Moh theorem can actually be assumed to be successive approximate 
roots of the equation F, and they can be computed by the algorithm of approx- 
imate roots. More precisely, one has 

,5,; = Ip{F, Fi) = degxResY{F, F,) ior 1 < i < h 



^Notice that ni<5i is always a multiple of So . Thus the three properties of the Abhyankar- 
Moh theorem are trivially satisfied if /i < 1 and d/i_|_i = 1. 



*In such semigroups, property (III) in proposition 3.6 is substituted by UiSi < (5i_|_i 



where Fi — app{di,F), Ip{F,Fi) denotes the intersection number at P of the 
projective completions of the curves given by F and Fi, and where Fi and Si 
are given by the algorithm of approximate roots (see H). In particular, one has 



—v-p{Fi) = 6i . Thus, algorithm 3.4 computes at the same time the generators 
provided by the Abhyankar-Moh theorem and functions achieving such values 
in an effective way. 

Even more, if the algorithm succeeds, i.e. if one arrives to the end with the 
properties required by proposition 3^ j, we are sure that there is one branch 



at P, and hence the curve is absolutely irreducible. In case of fail, i.e. if such 
conditions are not satisfied in any of the steps of the algorithm, one can conclude 
that the plane curve has more than one branch at P. 

Remark 3.9 In case of having a priori a parametrization for the singularity 
of the unique branch at P, we could use it to compute contact orders instead 
of using resultants because of the formula Si = Ip{F,Fi). This would give an 
alternative for the algorithm. 

Example 3.10 Consider the affine plane curve Y^ + Y^ +X'^ = defined over 
W2, with only one point at infinity P = (1 : : 0). The degree of the curve is 
multiple of the characteristic, so with the change X = X + Y'^, Y = Y one gets 
the plane model F(X, Y) = Y^ + Y^ + XY^ + X^Y^ + Y^ + X^, and one can 
apply the algorithm of approximate roots to F: 

Po=X,<5o=rfi=9,Fi=Y 

(5i = degxResY{F,Y) == 3 , da = gcd(9, 3) = 3 

F2 = app{3, F)^Y^+Y^+Y + X+1 

S2 = degxResY{F, F2) ^ 8 , ds ^ g cd {9,3,8) ^ 1 

Thus h = 2 and Sp = (9,3,8). As a consequence, there is only one branch at 



infinity since properties (I), (II) and (III) from proposition 3.t are satisfied. 



On the other hand, with the notations as in proposition 3.1, take a IF2 -&asis 
for A/B: 

Y(l+Y6) Y(1+Y6) 



hi 



X + Y3 " (X + Y3)(Y2+Y+1) 



X2+Y6 _ Y2(1 + Y3)(Y2+Y+1) 

h'\ — — iLa — r 

Y2 + Y + 1 X + Y3 

The values of this functions are —vp(hi) = 13 ^ Sp, — Wp(/i2) = 7 ^ T\,, 
—vp{h^) = 10 ^ Fp and —vp{h4) = 13 G Fp. Then change h^ by 

Y(1+Y3)(Y2+Y+1) 
g^^hi + hi^ X + Y3 

^Notice that the properties given by the Abhyankar-Moh theorem are just the same which 
are required by the criterion for one branch at infinity. 



10 



and now —vp{g4) = 10 G Fp. Thus, one still has to take the function 

Y(l + Y3)(Y4 + Y2 + 1) + (X + Y3)3 
04 = n.4 + n.1 + /icj = ■ — 

y4 4-r IT 3 (X + Y3)(Y2+Y+1) 

and now —vp{g4) = 4 ^ Fp. Hence, the Weierstrass semigroup at P is 
Tp = {0, 3, 4, 6, 7, 8, 9, 10, 11, 12, 13, 14, . . .} 

3.3 Application to codes 

In particular, the above results allow us to compute a basis of the vector space 
C{mP), for every m G Tp , by collecting just one function with an only pole 
at P of order r, for each < r < to with r G Fp . This is essential for the 
construction and decoding of one-point AG codes, and it can be easily done 
from the above results as follows. 



If r G Sp, one gets in an effective way the writing [•] from the remark 3.8 
in the form r ~ Ao(5o + -^i^i + . . • + Xh^h ■ Thus, if we have functions fi with 
iJp{fi) = ^^i (for instance, those which are obtained from the algorithm of 
approximate roots), then fr — fo° ' fi^ ' ■ ■ ■ ' fh'^ has an only pole at P o f ord er 
r. Otherwise, if r G Fp \ Sp the function fr is constructed by algorithm |3.2| . 

Remark 3.11 // we fix previously a semigroup with properties as in proposition 



3.6, we can try to find a plane curve F with an only branch at infinity achieving 
this semigroup, even if p divides gcd (^oi'^i)- 

This last procedure is just the inverse of that we have explained in this sec- 



tion, i.e. let the numbers Sq,6i, . . . ,Sh with the properties as in proposition 3.t 
be given, then one can construct in a recurrent way a sequence of polynomials 
which are called "approximants" of a polynomial F such that F = is a curve 
with a unique branch P at infinity and semigroup Sp generated by So, Si, ... ,dh 
(this is shown by Reguera in ^Ig] , where one can check more details). This has 
the advantage that one can built directly the examples together with functions 
(the approximants) having the generators of Sp as pole orders. 



Remark 3.12 By using the affine algebra B instead of the normal affine algebra 
A (with the notations as in section 3.1J, one can construct "improved Algebraic 
Geometry codes" as in UW. More precisely, such codes would be the dual of the 
functional codes given by lm{evpi{C{mP) fl B)). They can also be decoded by 
the Feng and Rao method in the same way. Now one only needs the semigroup 



Sp instead ofVp, so one avoids the triangulation algorithm 3.i. Moreover, 



by using remark 3.11 one can directly construct curves with only one branch at 



infinity and functions achieving the values of a semigroup Sp fixed a priori. 
Finally, the results of this section can be summarized as follows. 

11 



Theorem 3.13 Let T be an absolutely irreducible projective plane curve of de- 
gree m with only one branch P at infinity which is rational over the base field 
W. Assume the characteristic ofW is either or it does not divide simultane- 
ously m and m — ep(T). Then, by combining algorithms \3.^ and \3.^ one can 
compute the Weierstrass semigroup Tp and functions achieving the pole orders 
in Tp . As a byproduct, one gets a basis of the vector space CirP), for every 
r G T p . Moreover, for each semigroup and generators with the properties as 
in the Abhyankar-Moh theorem, one can generate curves with Sp equal to that 
semigroup, without restrictions on the characteristic. 

D 

In the next section, we will show how to calculate the Feng-Rao distance 
for Weierstrass semigroups which have been computed with the above method. 
In fact, we will focus on a general situation in arithmetic semigroups, just tak- 
ing into account the arithmetic properties of the semigroups in the theory of 



Abhyankar-Moh and the modifications given by proposition 3.1. 



4 Computing the Feng-Rao distance 

In this section, we will compute the Feng-Rao distance as a function defined on 
any arbitrary numerical semigroup, i.e. a subsemigroup of IN. Some formulae 
will stand for the general case, and then we will focus on some concrete types 
of semigroups, e.g. symmetric, Abhyankar-Moh (or, more generally, telescopic) 
and Weierstrass semigroups obtained as semigroups at infinity. Notice that such 
types of semigroups are also interesting in singularity theory. 

4.1 Apery systems and Feng- Rao distance 

In the sequel, we consider numerical semigroups, i.e. subsemigroups S* of IN 
such that tt(IN \ 5) < oo and G S". The number g = tt(]N \ S) is called the 
genus of the semigroup S. Since the genus is finite, there exists a (unique) 
element c £ S such that c — 1 ^ S and c + Z G S" for alH G IN. The number c 
is called the conductor of S, and one has c < 2g. Thus, the "last gap" of S is 
Ig = c — 1 < 2g — 1, where k is called a gap of S* if fc G INT \ 5. 

On the other hand, notice that every m > 2g is the (ttt, -f- 1 — g)-th element 
of S, that is pm+i-g , according to the notations of [|6|. Finally, the semigroup 

5 is called symmetric when r G 5 if and only if c — 1 — r ^ 5, for all r &1L. 
This is equivalent to say c = 2(7, that is to say Ig — 2g ~ 1. 

Definition 4.1 For any semigroup S* C IN with jJ(M \ 5*) < oo and G 5, the 
Feng-Rao distance of S is defined by the function 



JFR 



S — >TN 
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771 ^-^ dpRirn) = min{7^(r) \ r > m, r G S*} 

where v is the function 

V : S — yTN 

r ^ iy{r) = tl{(a, b) e S'^ \ a + b = r} 

With the above notations, the following result summarizes some known facts 
about the functions ly and Sfr. for an arbitrary semigroup. One can check the 
details in H or jl6|. 

Proposition 4.2 

(i) 6fr{0) = 1^(0) = 1, and 2 < i'{Tn),SFB.{m) < 777 + 1 tf m E S\{0}. 
(ii) ^{m) = 777 + 1 — 2g + D{m) for 777 > c rl where 

D{m) = tt{(a;, y) \ x,y are gaps of S and x + y — m} 

(iii) ly^m) — m + 1 — 2g for all m E S with m > Ag — 1. 

(iv) 6pF!(m) > 771 + 1 — 2g = d*{m — 1) for all m, G S, and equality holds if 
moreover m > Ag ^ 1. 

In particular, it follows that SpRim) = v{m) ~ m + 1 — 2g for all 777 G S' 
such that D{m) = 0. On the other hand, since from (iii) one has that v is 
an increasing function for m > 4g — 1 at most, the knowledge of the Feng-Rao 
distance is finitely determined by the values of ly. In fact, for any 777 G S* one 
has 

^FR.i'm) = min{7^(777), 7/(777 + 1), . . . , 1/(777')} 

where 777' is the least element in S with m' > m and such that 7/(777') = 
m' + 1 — 2g. Thus, elements m' € S with D{m') = are interesting for calcula- 
tions. Other interesting kind of elements in S are those satisfying the following 
equality: 

SpRim) = min{r G 5 | r > 777 + 1 - 2g} [@] 

In practice, such elements exist and have the advantage that the Feng-Rao 
distance for them is easy to compute. 

Our aim is to give some formulae which will allow us to compute the Feng- 
Rao distance for an arbitrary semigroup by means of an algorithm, provided 
a suitable set of generators is given. The main tool which will be used was 
introduced by Apery in Q in order to study semigroups of curve singularities, 
and it is nothing but the Apery systems of generators and their relations. 

^Notice that one usually assumes that m > 2g — 2 in coding theory, and hence m > c. 
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Definition 4.3 Let S* C IN 6e a semigroup with ^(IN \ S) < oo and e 5*; for 
e E S \ {0} define the Apery set of S related to e by 

{ao,ai,. .. ,ae_i} 

where ai = min{m €^ S \ m = i (mod e)} for < i < e — 1. 

Usually one takes e as the multiplicity of S, that is, e = eg = iiiin(S' \ {0}), 
but actually it is not necessary. On the other hand, notice that one could remove 
Oo = since it does not add any information about the semigroup. In the sequel, 
the index i will be identified to the corresponding element in Z/(e). In fact, 
one has a disjoint union 

e-l 

S=\J{a^ + eIN) 

i=0 

and therefore the set {ai, . . . , ag-i, e} is a generator system for the semigroup 
S, called the Apery (generator) system of S related to e. 

Moreover, let i,j G 2/(e) = Zg and consider z + j e Zg ; then 

Gi + Uj = ai^j + Oiije 

with Qfi.j > 0, by definition of the Apery set. The numbers aij are called Apery 
relations. 

Under these conditions, every m € S can be written in a unique way as 
m — Oi + le, with i G Zg and / > 0. Thus, we can associate to m two Apery 
coordinates (i, Z) G Zg x IN. 

In order to compute i'{m), set m = {i,l), a = {ii,li) and b = (12,12)', since 

m^ a + b = ajj+.i2 + {li + I2 + C(i^^i^)e 

then li + I2 = I — ctii,i2 ■ Write ii ~ k and i2 = i — k; ii I < ak,i-k the equality 
m = a + b is not possible, and so we are just interested in the case cuks-k < '• 
Thus, for < i < e — 1 and /i > define 

sf ) = fi{afc,»-fe <h\keT,} 

With these notations, the following result gives us a formula to compute ^{m). 
Proposition 4.4 v{m) = Bf^ + sf ^ + . . . + sf ' 

Proof : 

Suppose ak,i-k — h < I; then it has been considered at the right sum in the 
sets defining sf \ fif +^\ . . . , sf \ that is I - h + 1 times. 

On the other hand, the equality I1 + I2 = l~oik.i-k holds for l — h+\ possible 
pairs li,l2 , and so the theorem is proved. 

D 
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Now, if we want to have a formula to compute the Feng-Rao distance, the 
main remark is that i^(m) is increasing in I, because of the previous formula. 
Then it suffices to calculate a minimum in the coordinate i, what gives only a 
finite number of possibilities. More precisely, one obtains the following result. 

Theorem 4.5 With the above notations, set m = Oi + le. For each j G 
TLe , take ruj = Oj + tjc, where tj is the minimum integer such that tj > 



max < + ^, > . Then one has 

<^f_r("t-) — min{v{mj) \ j £ Zg} 
Proof : 

By using the above remark on the number v{m)^ the formula follows from 
the fact that ruj is the minimum element of S with first Apery coordinate equal 
to j such that ruj > m. 

a 

As a conclusion, computing the Feng- Rao distance is easy if we have the 
Apery set related to any non-zero element of the semigroup, and the method 
works in a quite general situation. Next, we will show how the above facts can 
be done more precise for the case of symmetric semigroups. 

4.2 Symmetric semigroups 

Now we will compute the value of SpRirn) for many elements m (^ S, S being 



a symmetric semigroup, and improve the computation in theorem 4.5 for them. 
The underlying idea is searching for the values m G S such that either Spnirn) = 
v{m) = d*{m — 1) = m + 1 — 2g or the formula [@] is satisfied; this formula 
will be called "minimum formula" in the sequel. Our results will partially cover 
lacks in the results given in jig] . 

First of all, symmetry can be easily expressed in terms of Apery sets. In fact, 
let e be any non-zero element of the semigroup S, and consider the Apery set 
{oo, oi, . . . , fle-i} related to e. For this Apery set, consider the index iV £ 2e 
such that 

ajy = max{ao, ai, . . . , Oe-i} 

One can easily check that the last gap of 5' is just Ig = cat — e, i.e. qn = c— 1 + e. 
Thus, it is easy to check that S is symmetric if an only if 

tti + aj^^i = a^ Vi G Zg 

In this case, notice that Ig — 2g — 1, and thus ajv — 2g — 1 + e. By using this 



fact and the formula of proposition 4.4, one obtains the following result, which 



provides us with a formula of 5pji{m) for a certain range of values in S. 
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Theorem 4.6 Let S be a symmetric semigroup. Then one has 

SpRifn) = v{m) =m — lg = m + l — 2g = e 
for all m = 2g — I + e with e G S \ {0}. 
Proof : 

Take an arbitrary non-zero e € S, and let N be the index in Zg such that 
otv is the maximum Apery element related to e. We will first prove the formula 
for all ?n G S" such that m = N {mod e), i.e. for m = on + ^e, with I > 0. 



By proposition 4.4, one has 



iy[m) — ti^ + B^ + . . . + B^ 

where B)^ = '\\{ak,N~k < /i | fc G 2e}. But because of the equality ak + ajsi-k = 
om for all fc G 2e , one has oik.N-k = for all fc G 2e , and thus Bj^ — e for 
all h > 0. Hence, i/(to) — (/ + 1) e = le + (a^v — Ig) — "m — Ig , and the 
formula is proved for v{m). The formula for 5pji{m) follows immediately from 
the definition of the Feng- Rao distance and proposition [4.2| (iv). 

In particular, for the fixed e the statement holds for I — 0, that is, for 
m = On ^ Ig + e. But e is an arbitrary non-zero element in S, and thus the 
theorem is proved. 

D 



Theorem 4.6 shows in particular that the minimum formula [@] holds for 



half of the elements in the interval [c — 1, 2c — 21. But, in fact, this conclusion 



of theorem 4.6 is also true for many other elements in this interval. In order to 
show it, we forget for a moment of Apery systems, and for m Cz I = [c, 2c — 2] 
define two integers n, q given respectively by the equalities 

m = n + c — 1 , n + q ^ c — 1 

Then one has m — 2c — 2 — q, and n E S iS q ^ S. Moreover, from proposition 



4.2 (ii) and the fact that S is symmetric, it follows that 

v{m) ^ n + v{q) 
where i'{q) = for q ^ S. Thus, in particular, if n G S" one recovers u{m) 



what gives a different proof of theorem 4.6 



Assume now that n ^ S and let n' be the least integer in S such that n' > n. 
Consider the element q' given hy n' + q' = c— 1 and let S = S{q) = q — q' = n' — n 
be the distance of q (resp. n) to a gap (resp. non-gap) of S. Notice that 
^il) ^ So — 1, Co being the multiplicity of S, since the interval [n,n') consists 
of gaps and there are at most Co — 1 consecutive gaps in S. 
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Theorem 4.7 Let S be a symmetric semigroup and m E I ~ [c, 2c — 2]. With 
notations as above, one has 

SpRim) = min{ri + ^{q), n + I + v{q — 1), .... n + 5 — 3 + v{q — (5 + 3), n'} = 

= min{i'(m), v{m + 1), . . . , i>(m + 5 — 3),n'} 

In particular, one has dFR.{m) < n' . 

Proof : 



Since n' £ S one has i>{c— 1 + n') = n'; thus, from 4.2 (iv) and the definition 
of the Feng-Rao distance one has 

SFj^{m) = niin{i^(?7i), i^lm + 1), . . . , i'{m + 5) — n'} 

But iy{m + 6-l) = n'-l + v{q' + 1) > n', since g' + 1 £ S*. Moreover, 
v{m + 5-2) = n' -2 + v{q' + 2) > n' , since g' + 2 e S* \ {0} and therefore 
v{q' + 2) > 2 (notice that q' can be negative but only in the case n' — c, and in 
this case q' = —1). This proves the theorem. 
D 



Remark 4.8 According to the above theorem, the Feng-Rao distance can be 



computed by using the formula given in proposition {.^ and taking the minimum 
of only S — 1 integers. In order to carry out this computation is better, in 
principle, to use the first equality in the formula of \4-% since integers with the 
size of q are smaller than those with the size of in. Next corollary complements 



the result of theorem ^.6 



Corollary 4.9 Let S be a symmetric semigroup and m E I . With notations as 
above, the minimum formula [@] holds for m if and only if one has 

m > m 

for every q in the interval [q' + 3, g] . 

In particular, this is satisfied for the cases n €z S, S — S{q) — I or 2, and 
when (5 = 3 but q is not an irreducible element of S (i.e. a non-zero element of 
S which is not a sum of two non-zero elements of S). 

Proof : 



It foUows from theorems 4.6 and 4.7. Notice that ^{q) > 3 if (5 = 3 and q is 
not irreducible. 
D 



Remark 4.10 The condition in the above corollary is easy to check, in practice, 
for many values of m as, for example, those with small values of q. 
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However, by using Apery systems relative to any non-zero element e of S 
(for instance e = eo) it is always possible to check it. In fact, u{q) can be 



computed from proposition ^.4 and, in order to compute S = 5{q), we can use 
the following result. 

Proposition 4.11 Assume that q G S, with Apery coordinates (related to e) 
given by q — ai + le. Then S is the least integer such that q — S < ai^g . 

Proof : 

On needs to find the value of S such that {q — S,q] C S and q~ 5 ^ S. But 
if g = i (mod e) one has q — j = i — j (mod e) for aU j G Zg . Thus, the 
proposition follows from the definition of Apery elements. 

D 

Finally, let qo be the smallest integer in S such that i^iqa) < S{qo), and set 



qo = c— 1 if such an integer does not exist. From corollary 4.9 one has go ^ eo + 2 
and qo ~ Co + 2 if and only if eo, gq + 1, cq + 2 G S and eg > 2. This happens, 
for instance, for the symmetric semigroup whose gaps are 1,2, . . . , g — l,2g — 1, 
with g > 3. For a general symmetric semigroup, the value of qo is, in practice, 
much bigger. 

Now, from proposition ^^ (iv), the minimum formula holds in the interval 
[Ag — 1, 00). In fact, from the above results one deduces the whole interval where 
this formula is satisfied. 

Corollary 4.12 For a symmetric semigroup S, the minimum formula [@] is 
satisfied in the interval (mo, 00), where tuq = 4g — 2 — qo , but it does not hold 
for mo . 

Proof : 



It follows from theorem 4.7, corollary 4.£ and the definition of qo ■ 
D 

Remark 4.13 Notice that mo > 4(7 — eo — 3, what gives an estimate for mo . 

Example 4.14 

• For the semigroup S = (9, 12, 15, 17, 20, 23, 25, 28) one has c = 32 and 
qo — 25, thus [@] is satisfied for m. > 38. 

• For the semigroup S = (6, 8, 10, 17, 19} one has c = 22 and qo — 19, thus 
[@] is satisfied for m > 24. 
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For the semigroup S = (8, 10, 12, 13) one has c = 28 and qo — 25, thus 
is satisfied for to > 31. 
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• For the semigroup S — (6, 10, 15) one has c = 30 and qq = 29, thus [@] is 
satisfied for m > 30. 

4.3 Semigroups at infinity 

In the situation of semigroups at infinity being computed as in section 3, we 
want to calculate the Apery set of the semigroup Sp related to e = ^o = deg T, 
with the assumptions of the Abhyankar-Moh theorem. In general, we can solve 
this problem for an arbitrary telescopic semigroup S. 

Recall that, for such a semigroup, the main fact is that any m (^ S can be 
written in a unique way in the form 



m = ^ XkSk = XqSq + ^ XkSk 

fc=0 fc=l 

with Aq > and < Xk < nk — dk/dk+i for 1 < k < h. But this means exactly 
that the possible Apery elements related to 6a are only those with Aq = 0. On 
the other hand, the number of all the possible elements of this form is 

nin2 ...nh = So 

since di — Sq and d^+i = 1. Thus, all of them are different modulo Sq and 
minimum with this property, and hence they are just the elements of the Apery 
set of 5* related to Sq ■ 

As a consequence, for a given m ^ S written in the form 



m = XqSq + ^ XkSk 

k=l 

we can compute its Apery coordinates I = Ao and i G 2/(5o) such that i = 
'^k=i'^k6k {mod do)- Thus, we can easily compute a^ and atj , and so iy{m) 
and SpRirn) for a telescopic semigroup and, in particular, for a semigroup at 
infinity Sp when computed by the algorithm of approximate roots, under the 
assumptions of the Abhyankar-Moh theorem. 

In fact, the values of 5f_r(to) are known for many values of m G S, i f S is 



telescopic. More precisely, apart from the results given by proposition O. if we 
assume moreover that Sh = maxjJo, Si, ... , Sh}, then the minimum formula [@] 
is true for every m G [c, 2c — 2] in the interval^ [(4g— 1) — {dh~l)5h,4:g — 1), i.e. 
one has qn > {dh ~ l)Sh . Also one has SpRim) = j + 1 if {j — l)Sh < rn < j dh < 
(dh — i)Sh (see Jlgl). Hence, the Feng-Rao distance is unknown in general for 
the interval [{dh — l)Sh + 1, (4^ — 2) — {d^ — l)Sh], but nevertheless it can be 

''^This result is improved by corollary |4.12| . For example, the telescopic semigroup S = 
{8, 10, 12, 13) the formula [@] is satisfied from m = 31 instead of m = 42, and for the telescopic 
semigroup S = (6, 10, 15) the formula [@] is satisfied from m = 30 instead of m = 44. 
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computed by the formulae of proposition 4.4 and theorem 4.5. Notice that, for 
example, one has 



J2^kSk] =(Ai-l)(A2-l)...(A,-l) 



\k=l 



for the Apery elements. 

Now, we want to compute the Apery set related to e = deg T for the Weier- 
strass semigroup Tp = Sp + &1IN + . . . + fe^IN, where Sp has been computed by 
means of the Abhyankar-Moh theorem, s = dim'p{A/B) and bi has been com- 



puted as in proposition 3.1. In order to do it, it suffices to solve the following 
question: 

For a given numerical semigroup S with Apery set {ap, ai, . . . , a^-i} 
related to e € S \ {0}, computing the Apery set {00,01, . . . ,0^37} 
for the semigroup S = S + 6IN related toe£5'\{0}. 

In fact, after having solved the above question, by repeating the procedure a 
finite number of steps one gets the Weierstrass semigroup Tp with its Apery 



basis (see remark 3.3) 



First, for a given b one wants to know wether b is already in S or not, because 
if it is so, then S ^ S + 61N and nothing changes. This can be easily tested with 
aid of the given Apery set of S as follows: 

(i) Calculate i G 1.^ such that i = b {mod e). 

(ii) Then, b e S ii and only if / = ^— ^ > 0. 

e 

On the other hand, the candidates to be the elements of the new Apery set 
are obviously the numbers TOj.a = aj + Xb with < j < e — 1 and < A < e — 1, 
i.e. the number of possibilities is small. Thus, we proceed as follows: 

1. Initialize 57 = o,; for 1 < i < e — 1 (obviously oq = oo =0). 

2. Take one of the elements mj,\ and compute its remainder i modulo e. 

3. Compare rrij^x with the value of a7; if 07 is greater we have to change its 

value, set ol = rrij^x and so on with the next mj,x ■ 

At the end of this procedure we get the Apery set for S, and the new Apery 
relations can be easily computed from them. 

Remark 4.15 The results of this section suggest to present the elements m in 
a semigroup S as couples (i,l) £ 2e x IN i/m = Oi + Z e. Thus, S is given by 
nothing but the data {e ; ai for < i < e — 1} with the restrictions aij > 0. 
In particular, the computations made in section 3.3 to find a function for every 
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element in the Weierstrass semigroup are easier if we use the Apery description 
of the semigroup and functions for each element of the Apery set. 

More precisely, one just saves the data {ao = 0, ai, . . . , ae-i, e} together 
with the associated functions {hg = 1, hi, . . . , h^^i, h^} satisfying —vp{hi) = ai 
for < i < e — 1, and —vp{he) = e. Thus, for any m Cz S we compute its 
coordinates {i, I) and then /„ = hi- hi, is the function associated to m. We save 
in general many generators, but the description of an element m (z S in terms 
of them is quite fast (just doing a sum and two divisions by e). However, if we 
save the sequence of non-gaps {pi} with the associated functions {fi} we may 
use much more space, since the ratio g/e can be as large as we want. 



On the other hand, the new elements added in the steps of the algorithm 3. 
modify the semigroup changing those data. Thus, such changes can be included 
in this algorithm, and so the semigroup Tp , the functions and the Feng-Rao 
distance can be computed simultaneously. In particular, this gives an effective 



answer to the idea which was proposed in remark 3.S. 

Finally, for the computation of the Feng-Rao distance we should take into 
account that in case of Tp being symmetri^ then the calculations can be sim- 



plified by means of theorem (.1 or its corollaries. 
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